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ABSTRACT 

Context. Relations linking the temporal or/and spectral properties of the prompt emission of gamma-ray bursts (here- 
C-h after GRBs) to the absolute luminosity are of great importance as they both constrain the radiation mechanisms and 

represent potential distance indicators. Here we discuss two such relations: the lag-luminosity relation and the newly 
^0 ' discovered duration-luminosity relation of GRB pulses. 

Aims. We aim to extend our previous work on the origin of spectral lags, using the duration-luminosity relation recently 
discovered by Hakkila et al. to connect lags and luminosity. We also present a way to test this relation which has 
originally been established with a limited sample of only 12 pulses. 

Methods. We relate lags to the spectral evolution and shape of the pulses with a linear expansion of the pulse properties 
around maximum. We then couple this first result to the duration-luminosity relation to obtain the lag - luminosity 
and lag - duration relations. We finally use a Monte-Carlo method to generate a population of synthetic GRB pulses 
HH ■ which is then used to check the validity of the duration-luminosity relation. 

Results. Our theoretical results for the lag and duration-luminosity relations are in good agreement with the data. 
They are rather insensitive to the assumptions regarding the burst spectral parameters. Our Monte Carlo analysis 
of a population of synthetic pulses confirms that the duration-luminosity relation must be satisfied to reproduce the 
O ' observational duration - peak flux diagram of BATSE GRB pulses. 

Conclusions. The newly discovered duration-luminosity relation offers the possibility to link all three quantities: lag, 
duration and luminosity of GRB pulses in a consistent way. Some evidence for its validity have been presented but its 
, origin is not easy to explain in the context of the internal shock model. 

Key words. Gamma rays bursts: general; Radiation mechanisms: non thermal 

1. Introduction and corrected for beaming in the Ghirlanda relation), the 

rr ^ 1 . . i-i ■ luminosity- variability relation (Reichart et al., 2001) illus- 

V> , The prompt emission of gamma-ray bursts is characterized trating the tendency of mminous bursts to be more variable 

2". . by the diversity of the observed temporal profiles. Some &nd the lag _ luminosity relation ( her eafter LLR) discovered 

G\ bursts show a simple shape with a fast rise followed by a by Norrig et &1 (20Q0) ^ Spectral , are & tQ quantify 

3 . slower decay while others have a complex structure with a the ch in the burgt profiles observed in differe nt energy 

succession of pulses which can be overlapping or separated bandg when yiewed &t u h pulseg arg narrower 

by intervals with almost no emission. Conversely the spec- and peak earUer Norrig et al (2Q00) cross _ correlated the 

tra are more uniform, generally wel fitted by two smoothly profiles between BATSE bands l and 3 and found that the 

connected power laws (the so-called Band spectrum; Band resulti k were decreasing with increasing burst peak 

et al, 1993) Many studies have tried to link the temporal luminosit As for the other re l at ions between luminosity 

£ . and spectral properties of bursts with the objective to gam and gpectral Qr temporal proper ties, the LLR offers clues 

insight into the physical processes governing the prompt tQ the ph ks of the prQmpt emission but also provides a 

emission. Already before and during the BATSE era sev- potential method to evaluate GRB distances from observa- 

eral relations between hardness and duration (Kouvehotou y ons at n i„ n enerK y only 

et al, 1993), intensity (Golenetskii et al, 1983) and flu- 8 8 
ence (Liang & Kargatis, 1996) were found. Following the 

discovery of the afterglows and the measure of the first Spectral lags are a direct consequence of the burst spec- 

redshifts, intrinsic quantities such as the luminosity or the tral evolution since a fixed, constant spectrum, would lead 

total radiated energy became accessible and new relations to proportional profiles in all energy bands. In a first paper 

appeared: the Amati (Amati et al, 2002) and Ghirlanda (Hafizi & Mochkovitch, 2007) we computed spectral lags of 

(Ghirlanda et al, 2004) relations between the peak energy pulses, defined as the time interval between pulse maximum 

of the global spectrum and the energy release in gamma- in two different bands. We obtained an explicit expression 

rays (assuming isotropic emission in the Amati relation for the lags and assuming the validity of an "Amati-like" 

relation between luminosity and the value of E p at pulse 
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Hakkila et al. (2008) have recently reconsidered the LLR 
and obtained a new relation which applies to individual 
pulses while the original Norris ct al. (2000) LLR consid- 
ered the burst as a whole. Moreover, they also found a 
correlation between pulse duration and luminosity (here- 
after DLR). These results offer the possibility to directly 
test and extend our previous work (Hafizi & Mochkovitch, 
2007). We start in Sect. 2 by comparing to observations our 
theoretical results for the LLR. Since they rely on the valid- 
ity of the DLR we propose to test it in Sect. 3 by comparing 
a synthetic population of GRB pulses to the observed peak 
flux - duration diagram of a sample of pulses collected by 
Hakkila & Cumbee (2009). We discuss our results in Sect. 4 
and Sect. 5 is the conclusion. 



2. The lag-luminosity relation 

2.1. Theoretical interpretation 

In paper I (Hafizi & Mochkovitch, 2007) we presented a 
simple analytical model to calculate spectral lags. We rec- 
ognized that lags were better defined using individual pulses 
rather than the whole burst profile. Pulses in the same 
burst can have different lags and Hakkila ct al. (2008) have 
shown that the global lag represents some average where 
the brightest pulse (which generally has the shortest lag) 
makes the dominant contribution. Looking for correlation 
between lag and luminosity it is therefore preferable to con- 
sider each pulse separately. Hakkila et al. (2008) obtained 



L = 6.110 52 (Aiia/O.OI s) 



-0.62 



erg.s 



(1) 



where L is the peak pulse luminosity and Aii 3 the spec- 
tral lag between BATSE bands 1 and 3. This expression 
is believed to be more directly linked to the physics of the 
prompt emission than the original one found by Norris et 
al. (2000) 



L = 1.3 10 M (Atis/0.01 s) 



-1.15 



erg.s 



(2) 



where lags were computed by cross-correlation of the full 
burst profile between the two bands. For individual pulses, 
spectral lags are more easily estimated from the time dif- 
ference between the peaks. These "pulse peak lags" gen- 
erally agree with those obtained by cross-correlation and 
have been used by Hakkila et al. (2008) to get Eq.(l). 

In our theoretical analysis (Hafizi & Mochkovitch, 2007) 
we calculated pulse peak lags from a linear expansion of the 
pulse shape and spectral properties around the maximum 
in BATSE band 1 (20 - 50 keV). Our result directly relates 
the lag to spectral evolution in a very transparent way. We 
get 



Aii 3 /i3,B p e p + /i3,a a + fi3,p b 



t„ 



Ci 



(3) 



where t p is the pulse duration. The other terms are defined 
in the following way: 



/13.X 



dLogFi?, 



dLogX 



with .F13 



(4) 



Here t\ is the time of pulse maximum in BATSE band 1 and 
Bapix) is the spectrum shape for which we assumed a Band 
function (with x = E/E p , E p being the peak energy of the 



spectrum and a and /3 the two spectral indices at low and 
high energy). The derivatives e p , d and b are respectively 
given by 



e =*t 



a 

a= — t F 

a 



(5) 



and are evaluated at t = t\. Finally C\ is a "curvature 
parameter" for the pulse around maximum. We have 



Ch _ iVifti) 



t 2 



(Ci < 0) 



(6) 



where N\(t) is the count rate in BATSE band 1. 

The functions fis,x depend on the spectrum param- 
eters E p , a and (3 at pulse maximum while e p , a and b 
represent their evolution. The two remaining quantities in 
Eq.(l), t p and C\, are fixed by the pulse shape. They show 
that spiky pulses (large |Ci|) have shorter lags than broad 
pulses (small |Ci|) for a given spectral evolution and pulse 
duration and that short pulses are expected to have short 
lags, both effects in agreement with observations (Norris & 
Bonnell, 2006; Gehrels et al., 2006; Hakkila et al., 2007). 
The three functions /i 3 .b p , /i3. Q and /i3,/3 have been rep- 
resented in Fig.l for a = —1 and (3 — —2.25 which are 
the central values of the distributions found by Preece et 
al. (2000) in their study of the spectral properties of bright 
BATSE bursts. Their behavior can be understood by noting 
that at large (resp. small) E p values T{E p ,a, p 1 ) depends 
on a (resp. (3) only. Therefore /i3,_e p = dLogTx^j 'dLogE p 
mostly contribute at intermediate E p (between BATSE 
bands 1 and 3) while /i3. Q (resp. f\s,p) dominates at large 
(resp. small) E p . The resulting ratio Atis/t p has been plot- 
ted in Fig. 2 as a function of E p for different values of e p , 
d, b and |d | = 10. 

This value of the curvature parameter as been adopted 
as representative of a "typical pulse" . In any case, the re- 
sults for a different |Ci| are easily obtained by rescaling 
Aii 3 /i p by a factor 10/|Ci|. As the maximum of E p gener- 
ally precedes that of the count rate in most pulses we have 
e p < 0. Similarly, the decrease of the spectral indices (spec- 
tral softening) begins before pulse maximum implying that 
d > and b > (since a and f3 are negative). 

To link spectral lags and luminosity Hafizi & 
Mochkovitch (2007) have moreover assumed an "Amati- 
like relation" between E p and the pulse peak luminosity of 
the form 



E n 



380 



1.6 10 52 erg.s - 



0.43 



keV. 



(7) 



This relation, proposed by Ghirlanda et al. (2005), is ex- 
pected to be valid at any time contrary to the original 
Amati relation (Amati et al., 2002) which applies to the 
burst as a whole. Using Eq.(3) and (7) it becomes possible 
to represent Atis/t p as a function of luminosity. The re- 
sults have also been plotted in Fig. 2. It can be seen that, as 
long as d and b are not too close to 0, Atis/t p only weakly 
depends on the luminosity. It is only if the spectral evolu- 
tion is limited to a decrease of the peak energy (the spectral 
indices staying constant) that Aii 3 /i p can reach very low 
values at both high and low luminosities. Since observed 
bursts show a simultaneous spectral evolution in E p and 
the spectral indices we expect that Aii3 / t p will not change 
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Fig. 1. Plot of the three functions /13 b p (full line), — /13 a (dotted line) and —fi3p (dashed line) for a = — 1 and 
£ = -2.25. 



much from pulse to pulse in agreement with observations 
(Hakkila et al., 2008). 

A roughly constant value of At\s/t p however raises a 
problem which was already mentioned in paper I. If pulses 
indeed satisfy a lag-luminosity relation with bright pulses 
having very short lags, some additional parameter has to be 
correlated to the luminosity. In paper I we tentatively pro- 
posed that pulse curvature could be such a parameter, lu- 
minous pulses being spikier and less luminous ones broader. 
But the recent discovery by Hakkila et al. (2008) of a possi- 
ble correlation between pulse duration and peak luminosity 
offers a new perspective which can naturally account for the 
LLR. 

2.2. The duration-luminosity relation of GRB pulses 

Hakkila et al. (2008) re-analized the seven BATSE bursts 
with known redshift considered by Norris et al. (2000) to 
establish the original LLR but they treated each pulse of 
these bursts separately. From the 12 selected pulses they 
obtained a new LLR (Eq.(l)) but also discovered an even 
tighter relation between duration and luminosity 

L = 3.4 1 52 t~ 85 erg.s" 1 . (8) 

Coupling this to Eq.(3) then provides a very simple way to 
get the LLR. An important difference with paper I is that 
we do not necessarily assume the validity of the "Amati-like 
relation" (Eq.(7)) to link E p to the luminosity. Even if E p 
and L are essentially independent quantities a LLR can still 
be obtained with the same value of the peak energy at all 
luminosities. In Fig. 3 and 4 we present the results for both 



the Ati3 - t p relation and the LLR. The full heavy line 
in each diagram corresponds to our reference case which 
assumes the validity of the Amati-like relation and uses the 
following values of the parameters: a = — 1, f3 = —2.25, 
e p = —0.5, a = b = 0.1 and \C\\ = 10. In addition we 
plot with thin lines a few other cases: (1) e p = —1 (other 

parameters unchanged); (2) e p = —0.25; (3) a = b = 
(other parameters unchanged); (4) a = b = 0.2; (5), (6) 
and (7) no Amati relation but constant E p = 250, 500 and 
1000 keV; (8) and (9) similar to 1) but with \C\ \ = 3 and 
30. All these lines define a narrow strip showing that both 
relations remain fairly robust even when the parameters are 
varied by large factors. 

Wc have also plotted in Fig. 3 and 4 the data points 
for the pulses belonging to the sample studied by Hakkila 
et al. (2008). It can be seen that the agreement with our 
theoretical results is satisfactory. 

3. A test for the duration-luminosity relation 

3.1. Method 

The results presented in the last section rely on the valid- 
ity of the duration-luminosity relation (DLR) for pulses. If 
confirmed, the DLR would also offer a new method to es- 
timate GRB distances, simpler and easier to use than the 
LLR (Hakkila, Fragile & Giblin, 2009). Bursts with sev- 
eral pulses give the possibility of multiple measures of the 
redshift, increasing the resulting accuracy. Conversely the 
identical redshift for all pulses in a given burst allows to 
test the DLR. Assuming a power-law of the form L oc t~ s 
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L (erg.s" 1 ) 
10 49 1 50 1 51 10 52 1 53 




E p (keV) 



Fig. 2. Ratio of spectral lag over pulse duration as a function of the peak energy at pulse maximum (lower scale) and of 
the peak luminosity (upper scale) assuming the validity of the Amati-like relation (Eq.(7)). The full lines correspond to 
e p = —0.5 and from top to bottom respectively to a — b = 0.2, 0.1 and 0. The dashed lines have the same a and b but 
e p = —0.25. The adopted value of the curvature parameter is \C\ \ = 10. 



Hakkila, Fragile & Giblin (2009) find s = 0.8 ± 0.4 for a 
sample of 53 multi-pulsed GRBs, which is consistent with 
the result obtained from bursts with known redshift. 

We have performed an alternative and independent test 
of the DLR using a synthetic population of GRB pulses 
for which we predict the resulting observational duration - 
peak photon flux (i° bs - P) diagram which is then compared 
to real data. The synthetic population is generated follow- 
ing a Monte-Carlo method similar to the one described in 
Daigne, Rossi & Mochkovitch (2006): for each pulse we 
draw a redshift z and a peak luminosity L. We then either 
link E p and t p to the luminosity with Eq.(7) and (8) or 
adopt log-normal distributions independent of L. We want 
to see if the predicted i° bs - P diagram is in better agree- 
ment with the data when the DLR is adopted. 

Even if t p and L are uncorrelated we expect a first trend 
purely due to cosmological effects as more distant pulses are 
globally weaker and have longer durations. If an intrinsic 
relation such as Eq.(8) is also satisfied the trend should 
be more pronounced. The analysis by Hakkila & Cumbee 
(2009) of a sample of pulses extracted from 106 long GRBs 
gives P oc tp ' 27 which is shallower than Eq.(8). The ob- 
served population is however affected by cosmological ef- 
fects (time dilation and k-correction) on the pulse duration 
(Norris, 2002) and by selection effects due to instrument 
threshold. We model these different effects to generate the 
simulated observational t° bs - P diagram from our synthetic 
GRB pulse population. 



3.2. The synthetic pulse population 

To get the distribution of the pulse parameters z, L, E p 
and t p in the synthetic sample we make the following as- 
sumptions: 

- Redshift z: as we only consider long GRBs which have 
massive progenitors, the GRB rate Rgrb(z) could a priori 
be expected to be proportional to the cosmic star forma- 
tion rate SFR(z). However recent studies (Daigne, Rossi 
& Mochkovitch, 2006; Guetta & Piran, 2007; Kistler et al., 
2008) have shown that at large z the GRB rate still in- 
creases while the SFR probably decreases or remains con- 
stant. This suggests that stellar populations at large z are 
more efficient in producing GRBs for reasons which are not 
well understood (reduced metallicity or/and IMF favoring 
massive stars). In this study we adopt a burst rate which 
follows SFR3 of Porciani & Madau (2001). This SFR which 
keeps rising at large z is not realistic as it would overpro- 
duce metals at early cosmic times but Rgrb oc SFR3(z) 
provides a good fit of the redshift distribution of Swift 
bursts. We then generate a table of N (N from 10 3 to 
10 6 ) values of the redshift and corresponding luminosity 
distance d\,(z) with z being distributed as 



dV SFR3(z) 

—L = TV d 



1+2 



d,Z ' f°° dv SFR3(z) 
JO dz 1+z 



(9) 



where is the comoving volume element in the concor- 
dance cosmology. 
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Fig. 3. Plot of the pulse duration as a function of spectral lag. The full thick line represents our reference model which 
adopts the Amati-like relation (Eq.(7)) and a = -1, /3 = -2.25, E = -0.5, d = b = 0.1 and \C\\ = 10. The thin lines 
correspond to nine other cases with different choices of the parameters (see text for details). The squares are the data 
points of the sample considered in Hakkila et al. (2008). 



- Luminosity L: we adopt a power law luminosity func- 
tion 3>(L) oc L~ s between L m ; n and L max - For the burst 
population it has been shown that a power law LF with 
1.5 < S < 2 can reproduce the Log N - Log P curve 
(Firmani et al., 2004; Daigne, Rossi & Mochkovitch, 2006). 
We adopt the same range of values here and vary L m j n and 
L max respectively from 10 50 to 10 51 erg.s -1 and from 10 53 
to 10 54 erg.s" 1 . 

- Spectral parameters: the peak energy is either ob- 
tained from the luminosity with the Amati-like relation 
or has a log-normal distribution of central value E p _ and 
width (T ■ When the Amati-like relation is adopted we add 
a dispersion a a around Eq.(7). We draw the spectral in- 
dices a and (3 in agreement with the distributions found by 
Preece et al. (2000) for bright BATSE bursts. In Daigne, 
Rossi & Mochkovitch (2006) the values of ctq, E Pt0 and a a 
were adjusted to provide a good fit of the E p distribution of 
bright BATSE bursts. We keep the same values as a starting 
point but we also vary them since we are now considering 
individual pulses rather than the entire bursts. 

- Duration: To get the pulse duration we either assume 
the validity of the DLR (Eq.(8)) with a dispersion cr tp 
or adopt a log-normal distribution of central value and 
dispersion adjusted to reproduce the observed distribution 
of pulse duration in the Hakkila & Cumbee (2009) sample. 

For a given pulse, the observable quantities P and i° bs 
are then computed from the intrinsic parameters. For the 



peak photon flux we have, assuming a normalized Band 
spectrum (i.e. J °° B(x) = 1) 

j / .300/_E° bs 

p = * 7 q i / B{x)dx (10) 

where E° hs = E p /(1 + z). To obtain the observed pulse 
duration we take into account both time dilation and the 
approximate dependence of pulse width with energy t p oc 
E~ 0A (Norris et al., 1996) so that 

t p bs ~(i + *)°-% . (ii) 

We adopt the threshold prescription from Band (2003) 
where the limiting photon flux is defined between 1 and 
1000 keV and depends on the observed peak energy. This 
finally allows us to construct the simulated observational 
i° bs - P diagram for comparison to real data. 

3.3. Results 

We define a reference case which corresponds to the follow- 
ing choice of the parameters: slope of the luminosity func- 
tion 5 = 1.7; L m in = 10 51 erg.s -1 ; peak energy obtained 
from the Amati-like relation (Eq.(7)) with an added dis- 
persion of 0.3 dex. We compare in Fig. 5 the resulting t° bs 
- P diagrams with and without the DLR. When the DLR 
is adopted we again assume a dispersion of 0.3 dex around 
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Fig. 4. Theoretical lag-luminosity relation for pulses compared to the data collected by Hakkila et al. (2008). The thick 
line is the reference case while the thin lines correspond to different model parameters (see Fig. 3 and text for details). 



Eq.(8). A fit of the diagrams by a power-law P oc t p s re- 
spectively gives s = 0.27 (with the DLR) and s — 0.09 
(without). In the first case, the agreement with the data of 
Hakkila & Cumbee (2009) is excellent while the correlation 
almost disappears in the second case. 

We then checked how these results are changed when 
we vary the model parameters and assumptions (see Table 
1): 

- Luminosity function: We list the power-law index s of 
the i° bs - P relation when we vary the lower and upper 
limits of the luminosity function i m i n and L max and its 
slope S. It can be seen that the results only weakly depend 
on L m in and i max and are nearly unsensitive to 8. 

- Peak energy distribution: we have first replaced the 
Amati relation by a log-normal distribution of central value 
E p = 600 keV and dispersion 0.3 dex which were the val- 
ues adopted by Daigne, Rossi & Mochkovitch (2006). Since 
these corresponded to the whole burst spectra and not to in- 
dividual pulses we have considered other possible E p values. 
We find that the power-law index of the £° bs - P relation is 
practically independent of the adopted E p , especially when 
we assume the validity of DLR. However in this case the 
t° bs - P relation becomes somewhat steeper than the data 
(with s = 0.33). 

- Dispersion of the DLR: we have increased the disper- 
sion of the DLR from er tp = 0.3 dex to 1.5 dex. We observe 
that the power-law index of the i° bs - P relation evolves 
from its reference value of 0.27 to 0.09 which corresponds 
to the situation without the DLR. It appears that the dis- 



persion of the DLR cannot exceed about 0.6 dex if we still 
want to fit the data. 



imin (erg.s L ) 


s (with DLR) 


s (without DLR) 


10 bi 


0.27 


0.091 


10 bU.b 


0.24 


0.093 


L 10bU 


0.22 


0.096 


imax (crg.S -1 ) 






1Q b3.b 


0.24 


0.093 


10 b4 


0.22 


0.096 


6 






2.0 


0.27 


0.081 


1.5 


0.27 


0.096 


E p (keV) 






400 


0.33 


0.089 


600 


0.33 


0.088 


800 


0.33 


0.096 


0t 






0.3 


0.27 




0.6 


0.22 




1.0 


0.15 




1.5 


0.09 




Threshold 






:2 


0.30 


0.11 


x2 


0.24 


0.085 



Table 1: Slope s of a power law fit (P oc t p s ) of the t p - 
P diagram with and without the assumption of the DLR 
(Eq.(8)) for pulses. In the six blocks we respectively vary 
the lower (i) and upper (ii) limits of the pulse luminosity 
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function; (iii) the slope of the luminosity function; (iv) 
the central value of a log-normal distribution for E p ; (v) 
the dispersion of the DLR; (vi) the detection threshold. 
The first row corresponds to our reference case with 
imin = 10 51 erg.s "\ i max = 1(F erg.s "\ S = 1.7, 
o"t = 0-3 dex. It also assumes the validity of the Amati-like 
relation (Eq.(7)) with a dispersion of 0.3 dex and adopts 
the threshold criterion for BATSE given by Band (2003). 
In each block only one parameter is varied, the others 
keeping the values corresponding to the reference case. 

Our study then indicates that an intrinsic correlation 
between pulse duration and luminosity is necessary to re- 
produce the observed i° bs - P diagram. This conclusion is 
not affected when we vary the model paramaters such as 
the luminosity function or the spectral properties of pulses. 
It is also robust regarding changes in the adopted sensitiv- 
ity: when we decrease (resp. increase) the threshold by a 
factor of two, the i° bs - P correlation becomes only slightly 
steeper (resp. shallower). 

Nevertheless, a few words of caution should be expressed 
since our analysis does not take into account possible ad- 
ditional selection effects which may not apply equally to 
pulses of different durations. For example the data has been 
collected with a trigger criterion applied to the full burst 
and not to individual pulses. It therefore includes some 
pulses below the threshold, coming from bursts which trig- 
gered at a brighter instant of the light curve. Also, the pulse 
selection and identification technique can fail when many 
pulses overlap, which is another source of selection effects, 
not easy to quantify. It is possible that these different biases 
contribute to produce an effective threshold with a limit in 
fluence in addition to the adopted limit in peak flux. A limit 
in fluence may artificially generate a trend in the P - f° bs 
diagram which could contribute to the observed relation. 

4. Discussion 

Assuming the validity of the duration-luminosity relation, 
we have tried to see if it can be understood in the con- 
text of the internal shock model for the prompt emission of 
GRBs (Rees & Meszaros, 1994). The isotropic luminosity 
generated by internal shocks can be approximated by 

L = Ef( K )e e (12) 

where E is the isotropic kinetic power in the relativistic 
flow, /(k) is the efficiency of dissipation by internal shocks 
(which mainly depends on the contrast k between the max- 
imum and minimum Lorentz factors) and e e is the fraction 
of the dissipated energy which is transferred to electrons 
and eventually radiated. 

It can be seen that the time scale r for variability of 
the Lorentz factor does not explicitely appear in Eq.(12). 
As the observed variability of the prompt emission reflects 
that of the Lorentz factor in the internal shock model (times 
the (1+z) dilation) any duration-luminosity relation im- 
plies that r should be in some way linked to E and pos- 
sibly also to k or e e . One could for example imagine that 
increasing E results in a more unstable outflow where the 
Lorentz factor fluctuates on a shorter time scale. This would 
induce a DLR which could become even more pronounced 
if the amplitude of the fluctuations (and therefore «) also 
increases with E. 



From an observational point of view there are some in- 
dications that E is anticorrelated with the opening angle of 
the relativistic jet (Frail et al., 2001). A thinner jet drilling 
its way through the envelope of the progenitor star could 
be more sensitive to Kelvin-Helmholtz instabilities devel- 
oping at its boundaries (Aloy et al., 2002). This would lead 
to a more irregular outflow with a shorter time scale of 
variability of the Lorentz factor, finally leading to a DLR. 
But clearly this discussion is somewhat speculative and it 
remains that the internal shock model does not provide by 
itself a simple and direct way to explain the DLR. 

Recently the internal shock model has also been criti- 
cized for a series of reasons such as its low efficiency, the 
difficulty to explain the standard value a ~ —1 of the 
low energy index of the spectrum (Ghiscllini et al., 2000) 
or the possible complete suppression of shocks if the flow 
is strongly magnetized. Proposed alternatives to internal 
shocks are reconnection processes (Giannos & Spruit, 2007) 
relativistic turbulence (Narayan & Kumar, 2009; Lazar, 
Nakar & Piran, 2009) or comptonized photospheric emis- 
sion (Bcloborodov, 2009). Unfortunately the modelling of 
these mechanisms has not reached a degree of accuracy 
where detailed predictions can be made on the properties 
of pulses. 

5. Conclusion 

We have considered the relations existing between spectral 
lags, duration and luminosity in GRB pulses. Extending 
a previous work by Hafizi & Mochkovitch (2007) we have 
first shown that the lag over pulse duration ratio does not 
vary much among pulses (remaining of the order of a few 
percents). This result holds as long as the spectral softening 
following pulse maximum is not limited to a decrease of 
the peak energy but also affects the spectral indices, as 
indicated by the observations. We have then included in our 
analysis the relation between pulse duration and luminosity 
recently discovered by Hakkila et al (2008). Combined to 
our results it allows to link all three quantities: At, t p and L. 
The lag-duration and lag-luminosity relations we obtain are 
in good agreement with the data. Also, they do not strongly 
depend on the assumptions for the spectral parameters at 
pulse maximum: values of E p , a, f3 and their derivatives, 
Amati relation or log- normal distribution of E p . 

Originally obtained with a limited set of only 12 pulses 
the DLR has recently received further support from the 
analysis of another sample of 12 pulses coming from 8 bursts 
detected by the HETE 2 satellite (Arimoto et al., 2010). Its 
validity however still needs to be confirmed and we have 
therefore proposed to test it in a different (statistical) way 
using the observational duration-peak photon flux (i° bs - 
P) diagram. For that purpose, we have adapted the Monte- 
Carlo code of Rossi, Daigne & Mochkovitch (2006) to gen- 
erate a sample of synthetic pulses for which we predict the 
observational i° bs - P diagram. It appears that the observed 
correlation P oc t p - 27 is reproduced only if pulses sat- 
isfy the DLR. This conclusion remains valid when we vary 
the pulse luminosity function and spectral properties (E p 
obtained from the Amati relation or having a log-normal 
distribution). Nevertheless we cannot completely exclude 
some bias in the pulse selection and characterization pro- 
cess which could contribute to the observed relation even 
in the absence of a DLR. 
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Fig. 5. i° bs - P diagram for a sample of 750 synthetic pulses corresponding to our reference case (see text for details). 
The left panel, which assumes the validity of the DLR looks very similar to the data collected by Hakkila & Cumbee 
(2009). The right panel, with a log-normal distribution of pulse duration do not show a significant correlation between 
the two quantities. 



We have finally confronted the DLR to the prediction 
of the internal shock model for the prompt emission. It ap- 
pears that the DLR cannot be obtained as a direct and 
simple consequence of the model. Additional assumptions 
are required, for example the possibility that the relativistic 
outflow becomes more unstable and variable when the in- 
jected kinetic power increases. Proposed alternatives to in- 
ternal shocks - reconnection, relativistic turbulence, comp- 
tonized photosphere - still don't have the predictive power 
to test if they can explain the DLR. 

In a future development of this work we plan to extend 
our analysis of the pulse properties (width and spectral 
lags) to other energy ranges. Data are sparse in the optical 
but the detailed light curve of the "naked-eye burst" GRB 
080319b has for example revealed interesting correlations 
between spectral lags at high and low energy (Stamatikos 
et al., 2009). At very high energy, Fermi observations have 
shown delays in the onset of the LAT component with re- 
spect to the MeV emission. Understanding the origin of 
these behaviors will provide clues for a better understand- 
ing of the prompt emission of GRBs. 
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